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ABSTRACT 



Aims. A high density portion of the Orion Molecular Cloud 1 (OMC-1) contains the prominent, warm Kleinmann-Low (KL) nebula 
plus a farther region in which intermediate to high mass stars are forming. Its outside is affected by ultraviolet radiation from the 
neighboring Orion Nebula Cluster and forms the archetypical photon-dominated region (PDR) with the prominent bar feature. Its 
nearness makes the OMC-1 core region a touchstone for research on the dense molecular interstellar medium and PDRs. 
Methods. Using the Atacama Pathfinder Experiment telescope (APEX), we have imaged the line emission from the multiple tran- 
sitions of several carbon monoxide (CO) isotopologues over the OMC-1 core region. Our observations employed the 2x7 pixel 
submillimeter CHAMP+ array to produce maps (~ 300" x 350") of 12 CO, I3 CO, and C ls O from mid-7 transitions (7 = 6 - 5 to 8 - 7). 
We also obtained the 13 CO and C ls O 7 = 3-2 images toward this region. 

Results. The 12 CO line emission shows a well-defined structure which is shaped and excited by a variety of phenomena, including the 
energetic photons from hot, massive stars in the nearby Orion Nebula's central Trapezium cluster, active high- and intermediate-mass 
star formation, and a past energetic event that excites the KL nebula. Radiative transfer modeling of the various isotopologic CO lines 
implies typical Hi densities in the OMC-1 core region of ~ 10 4 - 10 6 cirr 3 and generally elevated temperatures (~ 50 - 250 K). We 
estimate a warm gas mass in the OMC-1 core region of 86-285 M Q . 

Key words. Interstellar medium: nebulae, ISM: molecules, HII regions, Radio lines: ISM, Submillimeter: ISM, Radiative transfer 
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1. Introduction 

The Orion Molecular Cloud 1 (OMC-1) is a complex region of 
the interstellar medium (ISM) stretchin g over more t han 2.4 pc 
(20', roughly north-south) on the sky dKutner et al.L 1 1976b . Its 
densest part, toward which the Great Orion Nebula (M42), a 
classical compact "blister" Hn region, and its associated Orion 
Nebular Cluster (ONC) appear in projection, is one of the best- 
studied regions in astronomy. It is a test bed for studies of 
(proto)stars and clusters and the formation of low-, intermediate- 
, and high-mass stars (for a review of this region, see |0' Dell, 
■ I200U lO Dell et al.L 120081) . Much of this region's prominenc e 
is due to its distance of just 414 ± 7 pc dMenten et all 120071) . 
which makes the ONC and OMC-1 the closest regions of re- 
cent (few million years old) and ongoing high-mass star forma- 
tion. In the following, we shall use the term "OMC-1 core" or 
even just OMC-1 for the roughly 8' x 8'- or 1 pc 2 -sized dense 
molecular cloud region, which is located closely (0.1-0.2 pc) 
behind M42 and most of the stars in the ONC (see lZuckermanL 
119731: iGenzel & Stutzkil 119891 the latter give a comprehensive 
overview of this region and its phenomena). 

The OMC-1 core region may be divided into three main 
zones, all of which show bright (sub)millimeter wavelength 
emission from warm dust and molecular gas: the Becklin- 



Neugebauer/Kleinmann-Low (BN/KL) region, Orion South 
(OMC- 1 S or Orion-S) and the Orion Bar. 

The Orion BN/KL and Orion South regions are considered to 
be "twin" high-mass star forming regions because they have sim- 
ilar masses ($100 M Q ) and bolometri c lum inosities (10 4 - 10 5 



L Q , iMezger et all [1990; Dra patz et al. , 19831) and show co mpa- 
rable levels of star-forming activity (lO'Dell et all l2008[). The 
BN/K L region harbors the eponymous "hot core" (Masson et al., 
1984), which was (and commonly still is) taken to be the pro- 
totype of the hot dens e regions observed around many newly 
formed stars (see, e.g.. IKurtz et all [20001) . An interesting alter- 
native explanation for this region's energetics (other than being 
powered by an embedded central heating source) is a protostellar 
merger event th at released a few times 10 4 7 erg of energy about 



merger event mat released a tew tim es IIP erg of energy about 
500 years ago dBallv & Zinneckerl l2005t IZapata et all 120111: 
Ball v1it"ail 1201 lb - 
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The Orion Bar is a well-described photon-dominated region 
(PDR) located at the side of the OMC-1 core, facing M42, which 
is heated and partially ionized by far-ultraviolet (FUV) photons 
from the young massive stars (most of them from 6 l C, a spectral 
type 05-07 star ) that form the "Trapez i um" a t the center of the 
ONC (see, e.g., iHollenbach & Tielensl 119971: IWalmslev etail 
2000). In addition, the Orion Bar appears to be located at the 
edge of the Hn blister tangential to the line of sight. 

On giant molecular cloud (GMC) scales, low-rotational level 
(7) 12 CO emission (commonly from the 7=1-0 line) of the 
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Table 1. Observational parameters 
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Molecule/Line 


Frequency" 


E" / 

'-'Lip/ K 




^MB 


Receiver 


rW V 


r sys 






(MHz) 


(K) 


(") 






(mm) 


(K) 


c ls o; = 3 


-2 


329330.553 


32 


18.0 


0.75 


APEX-2a 


~ 1.0 


220 


"CO J = 3 


-2 


330587.965 


32 


18.0 


0.75 


APEX-2a 


~ 1.0 


210 


c 18 oy = 6 


-5 


658553.278 


111 


9.0 


0.47 


CHAMP+ 


< 0.7 


1600 


12 CO 7 = 6 


-5 


691473.076 


116 


8.6 


0.47 


CHAMP+ 


£0.5 


1800 


12 CO J = 7 


-6 


806651.806 


155 


7.4 


0.45 


CHAMP+ 


£0.5 


3900 


13 C07 = 8 


-7 


881272.808 


190 


6.7 


0.45 


CHAMP+ 


< 0.7 


3700 



Notes. Columns are, from left to right, CO isotopologue and transition, frequency, energy of the transition's upper energy level above the ground- 
state, the HPBW beam size, the main beam efficiency, the receiver used, the average precipitable water vapo r column and the a verage system 
temperature during the observations. (fl) Taken from the Cologne Database for Molecular Spectroscopy (CDMS, Miiller et al., 2005tQ- 



ambient, low density gas is usually used to trace the mass of the 
molecular ISM under a range of assumptions (Bloe men et al.L 
1984' lDame & Thaddeusl[T985h . However, in massive star form- 
ing regions with much higher densities and temperatures, obser- 
vations of the submillimeter and far-infrared (FIR) wavelength 
mid- or high-/ transitions of 12 CO and its 13 C and ls O isotopo- 
logues are required for determinations of the gas temperature 
and density, usually in conjunction with Large Velocity Gradient 
(LVG) radiative transfer and PDR modelings. 

Much of the 12 CO emission from the OMC-1 core has been 
proposed to arise from the neut ral and partially ionized b ack side 
of the Hn blister, i.e., the PDR dGenzel & Stutzkilll989l) . The pi- 
oneering submillimeter observation s of the 12 CO 7 = 7 - 6, 6-5 
and 13 CO 7 = 7 - 6 transitions by ISchmid-Burgk et al l (11989) 
and lGraf et al.l (1 19901) revealed high density gas (n > 10 4 cm' 3 ) 
and elevated temperatures (T > 50 K) in all parts of the OMC-1 
core region. Obser vations of the even m ore highly excited 12 CO 
7 = 9-8 line bv iMarrone et al.l (120041) with the Receiver Lab 
Telescope (84" resolution) showed that the hot and broad veloc- 
ity emission in the line profiles arises mainly from the BN/KL 
region, while much of the narrower (~ 3-6 km s _1 wide) emis- 
sion arises from the PDR. This was als o observed in the spectr a 
ofthe I2 C07 = 9-8 transition taken bv lKawamura et al.l(l2002l) . 
which indicated warm molecular gas from the extended "quies- 
cent ridge" region, north of Orion BN/KL. Further more, the ob- 
servati ons of 12 CO 7 = 7-6 and 7 = 4-3 made by Wilso n et al.l 
(2001) together with LVG modeling suggest that the broad line 
widths of these lines from the BN/KL region are probably due to 
shock heating, while most of the narrow lin e extended emission 
can be explained by PDR models. Recently. iFuruva & Shinnapal 
( 2009) presented the BN/KL region images in the 12 CO 7 = 7-6 
line. 

Clearly, multi-transition, high angular resolution imaging 
studies are needed to disentangle the relative contributions of 
shock and radiative excitation to the CO emission, and to iden- 
tify possible energy sources. All of the spectral line observations 
described above report either single pointing toward selected po- 
sition s in OMC- 1 or limited mapping of spe cial areas such as the 
BardLis et al.Ul998l) or the BN/KL region (IFuruva & Shinnagal 
2009 ). The only exception s are the 12 CO 7 = 7-6 map by 
Schmid-Burgk et al.l(U98 9) and the I2 CO 7 = 4- 3 and 7 = 7- 6 
maps bv lWilson et al.l (12001). However, the first of these, which 
covers the whole OMC-1 core, was taken with the poor resolu- 
tion of the Kuiper Airborne Observatory. 

Here we present the first high (better than 10") resolution 
large scale maps of the OMC-1 core using the 2 x7 pixel submil- 
limeter CHAMP + array receiver on the APEX telescope. Seven 
of the receiver units cover the 620 to 720 GHz frequency range 
and seven the 780 to 950 GHz range, allowing imaging of two 



lines simultaneously. As we shall see in $2] we imaged the 
OMC-1 core in 12 CO 7 = 6-5 plus 7 = 7 - 6 in one run 
and 13 CO 7 = 8-7 plus C I8 7 = 6 - 5 in another. The fast 
mapping afforded by the 14 unit instrument resulted in a supe- 
rior image consistency. The mid-7 CO data were complemented 
by the maps of the 13 CO and C 18 Q 7 = 3 - 2 lines. 



2. Observations 

The observations were performed in 2007 November and 2008 
July with the 12 meter APEX telescope on Llano de Chajnantor 
in Chile (iGiisten et al.L l20"06FL 

Large-scale maps of the I3 CO and C 18 7 = 3-2 lines were 
obtained in the On-T he-Fly (OTF) observat ions with the APEX- 
2a facility receiver (Risache r et al.L [2006). This double side- 
band (DSB) heterodyne receiver provided receiver DSB tem- 
peratures of 60 K and typical system noise from 100 to 200 K. 
Maps of size ~ 300" x 350" were centered on the Orion BN 
source (a,5j 2 ooo = 05 h 35 m 14 s .16, -05°22'2175). For all spec- 
tra, position switching to an emission-free reference position 
at offsets (-500", 0") was used. The Fast Fourier Transform 
Spectrometers (FFTS, iKlein et ail 120061) were operated with a 
1000 MHz bandwidth and a ~ 0.1 km s _1 spectral resolution. 

The submm observations of the high-7 CO lines were 
performed with the dual-color heterodyn e array CHAMP + 
ilKasemann et all 120061; IGiisten et al.Ll2008l) . operating 2 x 7 el- 
ements simultaneously in the 450 yum and 350 fim atmospheric 
windows. We observed two setups: in 2007 November the I2 CO 
7 = 6-5 and 7 = 7-6 lines were mapped simultaneously, 
while in the second coverage in 2008 July, the 13 CO 7 = 8-7 
and C 18 7 = 6 -5 transitions were observed in parallel. The 
MPIfR Array Correlator Spectrometer (MACS) was used as the 
backend during the former run, and the central pixels of each ar- 
ray were always connected to the FFT spectrometers mentioned 
above. Resolutions were ~ 0.42 km s for MACS and ~ 0.05 
km s" 1 for the FFTS. For the 2008 July observations, MACS 
had been replaced with a new FFTS array backend, operating 
for each pixel with a 2 x 1.5 GHz bandwidth and 8192 chan- 
nels per FFTS. We mapped the OMC-1 core region in a size of 
6' x 8' (0.72 x 0.96 pc 2 ) in OTF slews in the right ascension 
direction (each ~ 3' long), dumping data every 4". Subsequent 
scans were also spaced by 4" in declination so that the Nyquist 
sampled maps for each pixel were obtained. The angular resolu- 
tion varies with frequency between 6'.'7 (881 GHz) and 970 (658 



2 This publication is based on the data acquired with the Atacama 
Pathfinder Experiment (APEX) telescope. APEX is a collaboration 
between the Max-Planck-Institut fur Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. 
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Fig. 1. The OMC-1 core 12 CO 7 = 6-5 integrated intensity [-25, +30] km s image overlaid with the 7 = 7-6 integrated intensity 
[-25, +30] km s _1 contours running from 300 to 1000 K km s _! in steps of 100 K km s _1 , and the subsequent contours are plotted 
from 1200 to 7200 K km s 1 in steps of 600 K km s -1 . The 12 CO 7 = 6-5 (black) and 12 CO 7 = 7-6 (red) spectra are shown 
for Orion BN, Orion South, Orion Bar, and Orion East. Both images have been smoothed using Gaussian profiles (a width of 9'.'1 
for 12 CO 7 = 6-5 and 7'.' 8 for 12 CO 7 = 7-6) from their original Nyquist-sampled images. The black stars mark the positions 
of the five Trapezium stars (0 l Ori A, B, C, D, and E). The original beam sizes are shown at the bottom left of the image for 12 CO 
7 = 6-5 and 7 = 7 - 6 in red (9'.'6) and white (8'.'2), respectively. 



GHz). The same reference position at offsets (-500", 0") was 
used. The on-source integration time per dump and per pixel was 
1 second only. With the above over-sampling strategy, all pixels 
of CHAMP + covered a given grid position at least once. 

Data calibration was performed every ~10 minutes with cold 
and ambient temperature loads. The data were processed with 
the APEX online calibrator, assuming an image sideband sup- 
pression of 10 dB. The CHAMP + observing parameters, the tele- 
scope beam sizes, and typical values for the Precipitable Water 
Vapor (PWV) are summarized in Table Q] Typical single side- 
band (SSB) system temperatures were well below 2000 K and 
4000 K for the low- and high-frequency array of CHAMP + , re- 
spectively. 

Telescope pointing was established by 13 CO and C ls O 7 = 
3-2 line pointings every hour on nearby R Lup for the APEX- 
2a observations, resulting in an accuracy of ~ 2". In the shorter- 
wavelength submm windows, pointing is more difficult because 
of the much reduced system sensitivities and the paucity of suit- 
able pointing sources, strong enough for peak-up on their con- 
tinuum and/or line emission. During the 2007 CHAMP + obser- 
vations, o Cet was used as a line pointing source in the 12 CO 
7 = 6-5 mapping, but with an on-sky distance of ss 50° to 
the OMC-1 core, the absolute pointing on the latter was not ac- 
curately corrected. The wider IF processors and backends that 



became available in 2008 improved the situation significantly: 
simultaneous observations of the 12 CO 7 = 6-5 line and the 
CH 3 CN 7t = 36 9 - 35 9 multiplet (near 660.6 GHz) became 
possible in the low-frequency array of CHAMP + . We assume 
that the highly rotationally excited (lower energy level E/k-133 
K above the ground state) methyl cyanide line emission does 
peak on the Orion Hot Core proper, as defined by the cen- 
troid of 92 GHz CH3CN 7t = 5* — 4* emission determined 
with the BIM A interferom e ter, i. e., a,5 J2 ooo = 05 h 35 m 14 s .48, 
-05°22'3076 dWilner et al.L |l994|). Thereby, we established a 
pointing reference local to the OMC-1 core complex. While 
those small pointing maps toward the Orion Hot Core were 
made, the 12 CO 7 = 7-6 transition in the high-frequency ar- 
ray was observed in parallel, and provided the link to the 2007 
observations of the main CO isotopologues. The central pixels 
between the two arrays are co-aligned (better than one arcsec- 
ond). In practice, we shifted the 2007 maps such that the 12 CO 
7 = 7-6 line wing emission from -50 to -30 km s _1 match the 
pointed 2008 observations. The overall accuracy in the higher-7 
CO observations with CHAMP + is estimated to be < 4". 

All spectra were converted to the main beam brightness tem- 
perature unit, 7 M b = T*Jj] MB (77MB = B eS /F eS ), using a forward 
efficiency (F e ff) of 0.95 and beam coupling efficiencies (B e ff) as 
determined toward Jupiter (Table [TJ. The latter is motivated by 
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Fig. 2. The OMC-1 core 13 CO 7 = 8-7 integrated intensity [+5, +15] km s" 1 image overlaid with the C ls O 7 = 6-5 [+4, +12] 
km s _1 contours running from 15 to 65 K km s _1 in steps of 10 K km s _1 , and the subsequent contours are plotted from 80 to 200 
K km s _1 in steps of 30 K km s _1 . The 13 CO 7 = 8-7 (black) and C 18 7 = 6-5 (red) spectra are shown for Orion BN, Orion 
South, Orion Bar, and Orion East at the same positions as in FigureQ] Both images have been smoothed using Gaussian profiles (a 
width of 7'.'1 for 13 CO 7 = 8-7 and 9'.' 5 for C I8 7 = 6-5) from their original Nyquist-sampled images. The black stars mark 
the positions of the five Trapezium stars. The original beam sizes are shown at the bottom left of the image for 13 CO 7 = 8-7 and 
C ls O 7 = 6 - 5 in red (7'.'5) and white (1071), respectively. 



the average size of the emission in velocity space (Fig. |4j. All 
data were reduced using the standard procedures in the CLASS 
and GREG programs from the GILDAS packag^l 

3. Results and discussion 

3.1. Large scale maps - overall morphology 

Figure[T]shows « 6' x 8'-size maps of the OMC-1 core in the in- 
tensity of the 12 CO 7 = 6-5 and 7 = 7-6 lines integrated over 
the local standard of rest (LSR) velocity interval from -25 to 
+30 km s _1 . The emission distributions of both lines, clearly re- 
solved in our observations, resemble each other closely and are 
consistent with the JCMT /SCUBA 850 fjm image obtained by 
iJohnstone & Baiivl(ll999l) . where several regions are discernible. 
The most prominent of which are: (a) the emission maximum 
centered on the BN/KL region, (b) a secondary maximum co- 
inciding with Orion South (OMC-1S), (c) a region » 3' east of 
the Trapezium, which is known as Orion East (OMC-1 E; see 
below), and (d), the prominent straight Orion Bar PDR. The bar 
actually appears connected to the Orion South region, forming a 
mirrored letter L shape. Furthermore, there are two minor max- 
ima at 1 f8 N and 3. '3 NNE of the BN/KL peak position. 

Apart from their spatial distributions, the regions above can 
be distinguished by their centroid LSR velocities, Vlsr, and their 
line widths, AV. The Orion BN/KL region displays a spectacu- 
larly broad line up to 200 km s _1 wide, centered around 6 km s _1 
with a narrow (AV « 5 km s _1 ) feature at Vlsr ~ 8 km s _I 
superposed. Orion South is characterized by ~ 6 km s wide 
lines, with higher values at positions with outflow activities. The 
12 CO 7 = 6-5 and 7 = 7-6 lines observed toward the Orion 
Bar have FWHM widths of ~ 4 km s _1 and are centered around 



http://www.iram.fr/IRAMFR/GILDAS/ 



Vlsr = 10 km s~'. The Orion East region has prev i ously been 
mapped by iHoude etafl d2004l) and lHerrmann et alj d!997l) . but 
otherwise has received relatively little attention. The observed 
line widths and LSR velocities of Orion East, as well the high in- 
tensity of its CO emission (7mb Z 160 K) clearly demonstrate its 
PDR nature. A ll our measurements a re in good agreem e nt wit h 
the results of ISchmid-B urgk et al. (1989), iGraf etal] fl 990). 
IWilson et ail (1200 llhlKawamura et al l(l2002l) . and lMarrone et al.l 
(l2004h . The measurements of the CO isotopologues in the OMC- 
1 core are summarized in Table [2] 

Figure [2] shows the integrated intensity of 13 CO 7 = 8-7 
and C 18 7 = 6 -5 maps in a similar area as 12 CO but in a 
smaller range of radial velocities (from +5 to +15 km s _1 for 
13 CO and from +4 to +12 km s _1 for C I8 0). Both maps also 
show a good correspondence but with the integrated emission 
from C ls O 7 = 6-5 fainter than 13 CO 7 = 8 - 7. The two 
isotopologues have fainter and less broader lines compared with 
the I2 CO line emission which is expected for optically thick lines 
(see the spectra in Figs.[T]and|2]l. The Orion Bar is well resolved 
in both lines and also shows a sharp edge. In the mid-7 images 
(Figs. Q] and the strong peak in the center of the Bar is also 
seen in the recent Hersch el/SPIRE high-7 12 CO and 13 CO data 
obtained by Habart etafl (l2010h . In addition, the line widths of 
I3 CO 7 = 8-7 and C I8 7 = 6-5 toward all four main regions 
are narrow, even in the outflow zones, i.e., Orion BN/KL and 
South. 

The I3 CO 7 = 3 - 2 and C 18 7 = 3 -2 integrated inten- 
sity maps are shown in Figure [3] These images were made in a 
similar area as higher-7 I2 CO and in a range of radial velocities 
equal to I3 CO 7 = 8 -7 and C 18 7 = 6 -5. The C 18 7 = 3-2 
emission is only present where the 13 CO 7 = 3-2 emission is 
strong. The 13 CO and C 18 line profiles at this rotational transi- 
tion toward the Orion BN/KL and South regions are very broad 
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Fig.3. The OMC-1 core 13 CO 7 = 3-2 integrated intensity [+5, +15] km s -1 image overlaid with the C ls O 7 = 3-2 [+4, +12] 
km s contours running from 24 to 80 K km s _1 in steps of 8 K km s _1 . The 13 CO 7 = 3-2 (black) and C ls O 7 = 3-2 (red) 
spectra are shown for Orion BN, Orion South, Orion Bar, and Orion East at the same positions as in Fig.Q] The black stars mark 
the positions of the five Trapezium stars. The beam size of 21'.' 1 is shown at the bottom left of the image. 

Table 2. CO isotopologues measurements in the selected sources of the OMC-1 core 



Molecule/Line 




Orion BN 






Orion South 










T peak 


Vlsr 


AV 


^peak 


Vlsr 


AV 








(K) 


(km s" 1 ) 


(km s" 1 ) 


(K) 


(km s" 1 ) 


(km s" 1 ) 


11 CO J = 


8 


-7 


63.8+2.3 


10.3+0.3 


9.0+0.3 


41.0+3.6 


8.1+0.3 


4.9+0.3 


l2 CO J = 


7 


-6 


168.5+2.2 


9.8+1.1 


55.5+1.1 


165.6+2.0 


8.7+1.1 


6.9+1.1 


12 CO J = 


6 


-5 


164.6+1.6 


9.5+1.3 


53.3+1.3 


187.5+1.4 


9.1+1.3 


6.4+1.3 


c 18 oy = 


6 


-5 


19.0+1.3 


10.0+0.4 


5.8+0.4 


10.2+1.4 


7.6+0.4 


4.4+0.4 


"CO J = 


3 


-2 


68.7±0.5 


9.6+0.3 


5.6+0.3 


62.1 + 1.1 


7.3+0.3 


4.5+0.3 


c'«o; = 


3 


-2 


12.4+0.7 


10.2+0.3 


3.0+0.3 


19.8+1.6 


7.0+0.3 


3.3+0.3 










Orion Bar 






Orion East 




"CO J = 


8 


-7 


75.9+3.3 


10.7+0.3 


2.7+0.3 


55.7+2.5 


10.6+0.3 


2.3+0.3 


l2 CO J = 


7 


-6 


179.9±2.6 


10.5+1.1 


4.0+1.1 


156.2+3.1 


10.5+1.1 


4.0+1.1 


12 CO J = 


6 


-5 


180.9±0.9 


10.4+1.3 


3.8+1.3 


156.7+1.2 


10.4+1.3 


4.7+1.3 


c 18 o; = 


6 


-5 


23.9±2.2 


10.9+0.4 


1.8+0.4 


8.2+1.3 


10.4+0.4 


2.2+0.4 


"CO J = 


3 


-2 


90.1+1.1 


10.8+0.3 


2.4+0.3 


41.6+1.1 


10.6+0.3 


2.3+0.3 


c l8 o; = 


3 


-2 


13.4+1.0 


10.8+0.3 


1.7+0.3 


3.9+1.0 


10.4+0.3 


2.3+0.3 



Notes. Temperatures shown here are corrected to the main beam temperature unit, and line widths are measured in FWHM. 



and show well-defined wings. The southwestern part of the Bar 
is especially strong in this lower-7 transition for both isotopo- 
logues, which is also shown i n the Hersche l/SPIRE high-7 12 CO 
and I3 CO images reported by Haba rt et al.l (120101) . 

3.2. Velocity structure 

The 12 CO 7 = 6-5 and 7 = 7-6 velocity channel maps 
are shown in Figure |4] These maps reveal a well-resolved and 
complex structure of the OMC-1 core at ambient velocities as 
well as higher radial velocities. The outflows that are located 
mainly at Orion BN/KL, Orion South, and toward the north of 



Orion BN/K L are clearly seen in these images (see Zapat a et al.L 
120061 120091) . e.g., V LS r=0-6 km s 1 for the BN/KL region and 
Vlsr=12-14 km s~' for the Orion South region. Besides, no 
clear north-south velocity gradients across the OMC-1 core re- 
gion are observed. 

The compact and warm structures toward the BN/KL re- 
gion traced by the 12 CO 7 = 6-5 and 7 = 7-6 lines, 
which display a very broad range of velocities (up to about 
+ 100 km s~'), are part of the enigmatic molecular outflow that 
seems to be produced by a violent explos ion during the dis- 
ruption of a massive young stellar system ([Bally & Zinneckerl 
120051: IZapata et all 120091: iBallv et all 1201 II) . Our observations 
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200 100 -100 



R.A. offset (arcsec) 

Fig. 4. Low- velocity channel maps of the 12 CO 7 = 6-5 line with the 7 = 7-6 contours in the OMC-1 core. The contours are 
running from 50 to 200 K in steps of 25 K. Orion BN is located at the offset (0", 0"), and white stars mark the positions of the five 
Trapezium stars. 



reveal that some faint filamentary structures are very likely asso- 
ciated with the high-velocity 12 CO bullets (Peng et al, in prep.) 
reported bv lZapata etail (l2009t) . 

3.3. Line ratios 

In Figures EHTJ we show the line intensity ratio maps between 
the CO isotopologues and the different rotational transitions. 

From the image of the I2 CO to C I8 7 = 6 - 5 ratio (Fig. 
upper panel), we can see the spatial distribution of high optical 
depth regions at different velocities. The north-south dense ridge 
from Orion BN/KL to Orion South is clearly seen at 9 km s _1 . 
This high optical depth ridge is similar to the filamentary struc- 
tures of the NH 3 emission dWiseman & Hoi 1 19961 [l998), but 
has some morphological differences in the north of the Orion 
BN/KL region. Besides, the straight shape of the Orion Bar is 
pronounced at Vlsr=10-11 km s , where three high optical 
depth regions are seen at the two ends and center of the Bar 
at Vlsr=10 and 11 km s _1 , respectively. In the lower panel of 



Figure|5] the C 18 7 = 6- 5 to 7 = 3- 2 ratios around 10-1 1 
km s" 1 show a gradient in the Orion Bar which goes in the direc- 
tion from the Trapezium stars. It indicates that C 18 7 = 6 - 5 is 
strongly excited at the edge of the Bar by UV photons from the 
Trapezium stars. 

Figure [6] shows the ratio between 12 CO 7 = 7-6 and 
7 = 6 - 5 at three different radial velocities of 5, 10, and 15 
km s _1 . The three panels show clear variations across the inte- 
grated line emission. It is interesting to note that at the cloud 
velocity of 10 km s _1 , there are high 7-6/6-5 ratios located 
very close to the position of the Trapezium stars toward Orion 
BN/KL. These gradients are likely produced by those massive 
stars that heat the molecular cloud, causing the stronger 12 CO 
7 = 7-6 line emission compared with the 7 = 6-5 line. A 
gradient is also seen in the Orion Bar which goes in the perpen- 
dicular direction to the Trapezium stars (Fig. [6] at Vlsr = 10 
km s _1 ), and some patches with higher 12 CO 7 = 7-6 bright- 
ness temperatures are seen inside the Bar or behind the ioniza- 
tion front. Some horizontal and vertical strips, artifacts from the 
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Fig. 5. Upper panel is the 12 CO 7 = 6 - 5 to C ls O 7 = 6-5 ratio channel maps (#hpbw = 10") at different velocities from 6 to 1 1 
km s _1 in the OMC-1 core. The edge cut and empty area are due to different image sizes and data blanking (2 cr level of the C ls O 
7 = 6—5 temperature). Lower panel is the corresponding C ls O 7 = 6- 5 to 7 = 3- 2 ratio channel maps (#hpbw = 20"). The 
black arrow in the 10 km s _1 image indicates the ratio gradient across the Orion Bar (see text). The Orion Hot Core (HC) position 
is marked, and Orion BN is located at the offset (0", 0"). The positions of the five Trapezium stars are marked in each map. 



CO 7-6/6-5 




200 100 -100 



R.A. offset (arcsec) 

Fig. 6. Velocity channel maps (5, 10, and 15 km s~') of the 12 CO 
7 = 7- 6 to 7 = 6- 5 ratio (6» H pbw = 10") in the OMC-1 core. 
The 10 km s image is overlaid with the 12 CO 7 = 6-5 contours 
as shown in Fig. [1] The empty area are due to data blanking (3 
cr level of the I2 CO 7 = 7-6 temperature). The dashed-lines 
indicate that these ratio gradients are likely caused by outflows 
or filaments. Black triangles mark the Orion Hot Core position. 



OTF mapping, are also seen in the Orion Bar and South regions. 
These artifacts can affect the intensity ratio by ~ 14% given the 
calibration errors of 10% for both 12 CO 7 = 6 - 5 and 7 = 7 - 6. 



However, some structures seen in the CO 7 = 7- 6 to 7 = 6- 5 
ratio maps, e.g., the dashed-lines in Figure|6j are unlikely due to 
the calibration error or scanning strips during the observation, 
and seem to be footprints of filaments or outflows. Some simi- 
lar filament structure s have been noticed in the dust continuum 
emission observed bv lJohnstone & Ballvl (Il999h . 

In the ratio map of C lg O 7 = 3 - 2 to 13 CO 7 = 3-2 
shown in Figure [7] a clear elongated and dense ridge is seen in 
the north-south direction, where the 13 CO 7 = 3-2 intensity 
is weaker in the Orion BN/KL region than at Orion South and 
the north of BN/KL. Therefore, the lower CO column density 
derived from the lower C ls O 7 = 3 - 2 to 13 CO 7 = 3-2 ratio 
toward Orion BN/KL may be misleading. Even though the 13 CO 
7 = 3-2 line is optically thick and may be self-absorbed, this 
will result in a higher C ls O 7 = 3 - 2 to 13 C O 7 = 3-2 ratio 
instea d of a lower ratio in Orion BN/KL. As Goldsmith et al.l 
(1997) pointed out, the column density calculated using low-7 
C 18 lines is accurate except for the temperature > 150 K, which 
is the case in the OMC-1 core. Hence, higher-7 CO observations 
are critical to determine the CO column density in the OMC- 
1 core region. In addition, already shown in Figure [6] the low 
C ls O 7 = 3- 2 to 7 = 6- 5 ratios indicate that the C 18 7 = 
6-5 intensity is strongly enhanced toward the Orion BN/KL, 
South, and Bar regions, and is likely due to their PDR nature. 
Nevertheless, shocks/outflows may also play an important role in 
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Fig. 7. (a) The OMC-1 core C 18 7 = 3-2 to 13 C0 7 = 3-2 inte- 
grated intensity ratio map (6*hpbw = 20") overlaid with the 12 CO 
7 = 6-5 contours (grey) as shown in Fig. [1] The black contour 
represents the C 18 0/ 13 CO abundance ratio of 0.122 (60/490). (b) 
The OMC-1 core C 18 7 = 3 - 2 to 7 = 6- 5 integrated inten- 
sity ratio map (6>hpbw = 20") overlaid with the 12 CO 7 = 6-5 
contours as shown in Fig.Q] Black triangles mark the Orion Hot 
Core position. 



the Orion BN/KL and South regions, leading to an extra heating 
in the cloud. 



3.4. Excitation temperature and density estimates 

The 12 CO 7 = 6 - 5 and C ls O 7 = 6-5 data are used here to 
derive the gas excitation temperature and density of the OMC-1 
core. In the submillimeter regime where the Rayleigh-Jeans ap- 
proximation is often not valid, the observed radiation tempera- 
ture in local thermodynamic equilibrium (LTE) can be expressed 



rr,, hV 



i 



i 



,hv/kT a _ l e hv/kT bs _ j 



(1 



T "). 



(1) 



where 7 ex is the excitation temperature and t v is the optical depth 
at a specific molecular line transition. The background bright- 
ness temperature 7b g includes the cosmic background radiation 
of 2.73 K and the radiation from warm dust, ranging from about 
5 K for less den s e gas to about 26 K for the Orion Hor Core 
dGold smith et al., 1997). The second term inside the brackets 
can be neglected since it only contributes < 2% of 7^. By mak- 
ing the assumption of the same excitation temperature for 12 CO 
and C 18 7 = 6 - 5, the optical depths of 12 CO and C ls O 
7 = 6-5 can be determined from the relation 



7*( 12 CO) l_ e ^( ,2 co) 
7*(C I8 0) ~ 1 - e -T(C s O) ' 



(2) 



Since the optical depth of 12 CO 7 = 6 — 5 » 1, the optical 
depth of C ls O 7 = 6 - 5 can be directly obtained (~ 0.08) by 
assuming that the optical depth ratio is approximated to the iso- 
topolog ic abundance rat i o, and here we adopted f 12 C O|/[C 18 Q1 
of 490 dBoreiko & Betzt Il996t IWilson & MatteucciL 1 19921) and 



a beam filling factor of unity (TZ = 7mb)- The excitation tem- 
peratures of 12 CO and C ls O 7 = 6-5 can be estimated via 



hv 

' ex — ~T~ 

k 
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hv 



fcr MB ( 12 co) 



Then the total C 18 column density can be derived from 



JV(C 18 0) = 



8tt 3 v5^ 2 



TymdV 



1.33 xlO 12 (7 ex + 0.88)e* 



T M B<iV cm 



-2 



(3) 



(4) 



(5) 



where E up is 



110.6 K for the C I8 7 = 6 - 5 transition and 
<2rot is the rotational partition function. In the end, the H2 col- 
umn density can be estimated by adopting a [ 12 CO]/[C 18 0] 
abundance ration of 490 a nd a 12 CO abundance of 8 x 10" 5 
(Wils on & Matteuccil[l992l) . 

The results of the H2 column density and excitation tem- 
perature distributions are shown in Figure [8] where the north- 
south dense ridge near the Trapezium cluster and the Orion 
Bar in the southeast are clearly seen. The average t(C 18 0) at 
the 7 = 6-5 transition over the whole OMC-1 core region 
is about 0.08 with an average excitation temperature of about 
115 K, wh ich are consisten t with the 12 CO 7 = 9-8 obser- 
vations by Kawam ura et alJ (120021) . Besides, the minimum ex- 
citation temperature in the OMC-1 core is ~ 30 K, indicating 
a generally warm environment. The average C 18 column den- 
sity in the OMC-1 core is about 8.7 x 10 15 cm" 2 , corresponding 
to an H2 column densit y of 5.5 x 10 22 cm" 2 assum ing a I2 CO 
abundance of 8 x 10" 5 (IWilson & Matteuccil \l99% . Therefore, 
the dense gas mass is ~ 140 M Q in an area of 0.13 pc 2 in the 
OMC-1 core. The H2 column density agrees with the density 
of 5 x 10 21 - 7 x 10 22 cm " 2 derived from the m id- 7 CO iso- 
topologue observations by IWirstrom et al.l d2006l). O ur results 

1 et al.1 1 



agree with the recent measurement by IWilson et al.1 d201 ll) us- 
ing the 13 CO 7 = 6- 5 line with a l ower r 12 CQ1/fH 2 1 ratio of 
2 x 10" 5 . Additionally. Hab art et all d2010l) derived an H 2 col- 
umn density of about 9 x 10 22 cm" 2 toward the Orion Bar from 
the Herschel/SPTRE high-7 CO observations, assuming a lower 
excitation temperature of 85 K and a higher 12 CO to C 18 ratio 
of 560. Their result is consistent with ours, where we obtained 
an H2 column density of as 10 23 cm" 2 toward the Orion Bar with 
a higher excitation temperature (> 150 K). 

The JCMT 870 nm dust continuum emission image 
(iJohnstone & BallvL 1 19991) shown in Figure [8] agrees with our 
density map well, where the dust emission traces mostly the 
dense ridge (Orion BN/KL and South) and Bar regions. It is 
interesting to note that the dense ridge has an offset from the 
peak temperature positions, especially in the Orion South region 
(~ 30") and the north of the Orion BN/KL region (~ 20"). This 
dense ridge seems to extend farther to the north of the OMC- 1 
core, an d is probably re l ated to the N H3 filamentary structure 
seen by [Wiseman & Hoi d!998l 119961) . Additionally, as Figure 
|8] (b) demonstrates, the interface between the ionized and the 
molecular warm dense gas is very pronounced, where Orion 
BN/KL and Orion South are located at the north-south dense 
ridge, and the Orion Bar and Orion East are part of the high 
temperature enclosed structure well illustrated by the free-free 
continuum emission. Furthermore, our large-scale images are 
consistent with the [C n] and [Oi] emission maps obtained by 
iHerrmann et al.1 (119971) . where the rather uniform [Cn] emission 
indicates that PDRs are present over the whole OMC-1 core re- 
gion. 
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Fig. 8. (a) The C 18 column density image of the OMC-1 core overlaid with black contours of the JCMT 850 /urn dust continuum 
image. The 850 jum image is taken from the JCMT data archive. The contours represent 0.2%, 0.5%, 1%, 2%, 4%, 7%, 20%, 40%, 
60%, and 80% of the peak intensity (87 Jy beam" 1 ), (b) T he excitation temperature image for 12 CO and C ls O overlaid with the 35 
mm VLA-GBT continuum emission (Dick eTet all 12009) in black contours running from 15% to 95% in steps of 8% of the peak 
intensity (0.94 Jy beam -1 ). The stars mark the positions of the five Trapezium stars, and black triangles represent the position of the 
Orion Hot Core. 



3.5. RADEX modeling 

The non-LTE radiat ive transfer program RADEX 
(Ivan der Tak et al.L 120071) was used to investigate the accu- 
racy of the temperature and density calculation in LTE shown 
above. The optical depth effects are taken into account by 
RADEX using the escape probability approximation, where a 
uniform sphere geometry was chosen. The different transitions 
of 12 CO, 13 CO, and C ls O were used in the modeling, and their 
line widths were fixed at 5, 5, and 4 km s _1 , respectively, which 
are the average values from our spectra. The molecular data 
using in the modeling are taken from the Leiden A tomic and 
Molecular DatabaseQ (LAMDA; ISchoier et all [2005). In addi- 
tion, only H.2 was chosen as a collision partner, and 2.73 K was 
adopted as the background temperature. The input parameters 
are kinetic temperature and H2 number density. Each model was 
iterated with a kinetic temperature ranging from 50 to 450 K 
and an H2 number density ranging from 10 4 to 10 8 cm . The 
input column densities of these three molecules were chosen, 
so that their ratios were fixed at their isotopologic abundance 
ratios, i.e., 490 for [ 12 CO]/[C 18 0] and 60 for [ 12 CO]/[ 13 CO]. 
Three models were iterated with three different 12 CO column 
densities of 1 x 10 18 , 5 x 10 18 , and 1 x 10 19 cirT 2 which 
were obtained from the LTE calculation. Since these models 
did not take outflows into account, the output molecular line 
radiation temperatures (7r) were directly compared with the 
peak temperatures instead of the integrated intensities to avoid 



4 http://www.strw.leidenuniv.nl/~moldata/ 



the strong line wing emission. The modeling results are shown 
in Figure [9] 

Figure [9] reveals that most regions in the OMC-1 core have 
an H2 number density of ~ 10 4 — 10 6 cm 3 . The Orion Bar 
and Orion East emission peaks (Fig.|9]a), however, show higher 
densities. The distribution of the H2 number density in the 
Bar is c onsistent with 10 4 - 1 5 cm" 3 for a homogeneous 
medium (IWvro wski et al], 1 19971) . The modeling results in the 
Orion BN/KL, South, and Bar regions all show similar trends of 
a diverse kinetic temperature (about 100-200 K), which are con- 
sistent with the LTE calculation results. In addition, the non-LTE 
modeling results indicate that the central region of Orion BN/KL 
close to the Hot Core (in 20" size) has a number density of 
about 10 5 cm" 3 and a temperature of about 200 K, which is also 
close to the excitation temperature calculated in the LTE case. 
Noteworthily, the Orion BN/KL region and part of the Orion 
South region (Fig. [9] a) show very high kinetic temperatures 
(> 350 K) which may indicate an extra heating mechanism, e.g., 
outflow/shock heating. The suggestion of outflow/shock heating 
is supported by the larger I2 CO line widths (Fig.Q~|l and a great 
amou nt of known outflows in these regions (see iHennev et all 
2007). However, the high kinetic temperatures (> 250 K) in the 
Orion Bar and Orion East regions (Fig.|9]a) are hardly explained 
by the outflow/shock heating for the lack of outflow activities 
or broad 12 CO line widths. Instead, as PDRs, the Orion Bar and 
Orion East are expected to be heated by FUV photons from the 
Trapezium OB stars. M oreover, a recent study of the Orion Bar 
([Pelle grini et all 120091) suggests an extra heating by the excess 
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Fig. 9. The RADEX modeling results shown in the 12 CO 7 = 6-5 peak temperatures and the I3 CO/C ls O ratios for three different 
12 CO column densities. The column densities of 13 CO and C I8 were fixed at the isotopologic abundance ratio of 1/60 and 1/490 of 
the 12 CO column density. The line widths of 12 CO, 13 CO, and C ls O were fixed at 5, 5, and 4 km s _1 , respectively. The red and blue 
contours denote the temperature and H2 number density, respectively, and the color points represent the data in the different regions 
of the OMC-1 core. Each point represents a peak 12 CO 7 = 6-5 temperature and a temperature ratio between 13 CO 7 = 8-7 and 
C I8 7 = 6-5 from a single pixel with a resolution of 20". (a) The modeling input column density of I2 CO is 1 x 10 18 cirT 2 , and 
the data are plotted in a I2 CO column density density of 7.5 X 10 17 - 2.5 x 10 18 cm -2 , (b) The modeling input column density of 
12 CO is 5 x 10 18 crrT 2 , and the data are plotted in a density range of 2.5 x 10 18 - 7.5 x 10 18 cirT 2 . (c) The modeling input column 
density of 12 CO is 1 x 10 19 cirT 2 , and the data are plotted in the density of > 7.5 x 10 18 crrT 2 . Upper right panel shows the 12 CO 
column density map toward Orion BN/KL in contours of 5, 10, 20, 30, and 40 x 10 18 crrT 2 . The stars mark the positions of the five 
Trapezium stars, and the Orion Hot Core (HC) position is also marked. 



density of cosmic rays, which are trapped in the compressed 
magnetic field. 

3.6. Mass of the OMC-1 core 

By adopting a radius of 40"-60" and an Eb density of 5 x 10 4 - 
2 x 10 5 cirT 3 for Orion BN/KL, we can estimate a gas mass of 
7-85 M with the assumption of a spherical volume at a distance 



of 414 pc dMenten et al.L 120071) . In a similar way, the gas mass 
in Orion South is 3-49 M Q with a 30"-50" radius and an H 2 
density of 5 x 10 4 - 2 X 10 5 crrT 3 . For Orion East, the H2 density 
is ~ 2 x 10 5 cirT 3 with a 15"-25" radius, and the gas mass is 2- 
7 M . In contrast, it is very difficult to estimate the gas mass in 
the Orion Bar because of its geometry. We adopted a cylindrical 
geometry and an H2 density of 5 x 10 4 - 2 x 10 5 cirT 3 as the 
lower limit since an edge-on plane may contain much more gas. 
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With a length of 300" and a radius of 20"-30", the gas mass 
of the Orion Bar is 9-79 M . The estimated total dense H2 gas 
mass in the OMC-1 core is 21-220 M Q within a radius of 0.3 pc. 
This is also consistent with the mass estimated (~ 140 M ) from 
the LTE calculation above. If we include the warm diffuse gas 
(10 4 cm" 3 ) in the same volume, the total warm gas mass is 86- 
285 M in the OMC-1 core. This mass estimate for the OMC- 
1 core is close to that derived from the Odin 12 CO and 13 CO 
/ = 5-4 observations dWirstrom et al.L 120061) . i.e., 320 M Q in 
the molecula r ridge for a farther distance of 500 pc. The result of 
IWilson et al.l d2001l) also suggests a similar estimate of a warm 
gas mass of 31 0^30 M from the observations of 12 CO/ = 7-6 
and 7 = 4-3, again for a distance of 500 pc. Scaling these two 
literature values to the distance of 414 pc, results in 213-295 M 
for the total mass of the OMC-1 core. This is comparable to the 
value derived by us. Hence, we conclude that the higher-/ CO 
lines do trace most of the gas in the OMC-1 core which is heated 
mainly by the Trapezium stars. 

4. Summary 

In this paper, we present the first large-scale images of highly ex- 
cited 12 CO, 13 CO, and C ls O lines in the OMC-1 core. The high 
excitation temperatures (~ 150-200 K) reveal a hot enclosed 
structure which is mostly heated by the Trapezium cluster. The 
LTE approximation and non-LTE RADEX program were used to 
derive the excitation temperature and H2 density in the OMC- 1 
core. We found a typical density of 10 4 - 10 6 cm 4 in this re- 
gion. The clear dense ridge in the north-south direction is seen 
and is offset from the high temperature enclosed structure. Orion 
BN/KL and Orion South are located at this dense ridge, and the 
Orion Bar and East are part of the high temperature enclosed 
structure. The estimated mass of total warm gas is 86-285 M , 
assuming different geometries for different regions. In addition, 
the higher-/ CO lines trace most of the molecular gas in the 
OMC-1 core. The detailed investigation of individual region of 
the Orion BN/KL, Orion South, Orion Bar, and Orion East will 
be presented in the following papers. 
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